Electric field dependence of spin coherence in (001) GaAs/AlGaAs quantum wells 
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Conduction electron spin lifetimes (Ti) and spin coherence times {T2) are strongly modified in 
semiconductor quantum wells by electric fields. Quantitative calculations in GaAs/AlGaAs quantum 
wells at room temperature show roughly a factor of four enhancement in the spin lifetimes at optimal 
values of the electric fields. The much smaller enhancement compared to previous calculations is 
due to overestimates of the zero-field spin lifetime and the importance of nonlinear effects. 
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PACS numbers: 72.25.Rb,72.25.Dc,85.75.Hh 

The emerging field of semiconductor spintronics con- 
cerns the encoding, manipulation and detection of co- 
herence in the spin degree of freedom of mobile electrons 
in semiconductors^ '^. Configuring semiconductor materi- 
als in quantum wells (QWs) or other artificial structures 
can shorten the spin lifetimes and coherence times^ or 
lengthen them4. Electron spin coherence times are sen- 
sitive to carrier densitiea^iSii and electric fields^. A se- 
ries of calculations has suggested that an electric field 
{J^) can dramatically lengthen spin lifetimes in (001) 
QWs for spins oriented along either the (110) or (110) 
directions2iili. Similar considerations have led to the 
proposal of coherent diffusive spin rotation in QWs^^'^^. 
These calculations rely on perturbative models of the spin 
splitting in zincblende semiconductor QWs, where the 
= splittings are the sum of linear and cubic terms in 
the electron crystal momentum K and the splitting pro- 
portional to the electric field (Rashba Hamiltonian)iii4 
is also linear in K. 

Here we report detailed calculations of spin life- 
times (Ti) and spin coherence times (T2) in several 
GaAs/AlGaAs QWs at room temperature as a function 
of electric field. These calculations are performed us- 
ing a fourteen-band electronic structure theorjii^ii^ which 
treats the spin splitting non-perturbatively (to all orders 
in K). We find that the Rashba spin splitting deviates 
from a linear dependence on K for energies within about 
100 mcV of the band edge. We further find that the 
room-temperature electric-field-induced enhancement in 
the spin lifetime, for QWs between 50A and 150A, is only 
a factor of four. 

The Hamiltonian for bulk semiconductors in the pres- 
ence of an external electric field is 
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where mo is the free electron mass, c is the velocity of 
light, e is the elementry charge, p is the electron momen- 
tum operator, V^(f) is the periodic crystal potential, (t 
is the Pauli spin operator, and Kxt(f) = • f is the 
scalar potential generated by the external electric field 



The Schrodinger equation for heterostructure super- 
lattices (SLs) is written as 



ff"^(r|/:,5,K) = Ecsm{v\C,S,K.), 



and the SL Hamiltonian H^^ is given by 
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where Hi is the crystal Hamiltonian of the ith layer 
[Eq. iQJ] with parameters tabulated in Ref. 0, iViaycr 
is the number of layers in the SL unit cell, |£, 5, K) is 
the SL eigenstate for a carrier with wavevector K, pseudo 
spin quantum number S G {t, i}, and band index C, and 
i?£5(K) is the corresponding SL eigenenergy. To ob- 
tain \C, S, K) and EcsC^), we first solve the Schrodinger 
equation [Eq. at the zone center in which the SL 
wavevector K = 0. The solution, 



(r|7V,5,0) =^ ^ArSn(T(r)(r|71, cr, 0), 
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where |n, cr, 0) are the corresponding zone-center Bloch 
states of the constituent bulk semiconductors and the ex- 
pansion coefficients F^sna (r) are slowly varying envelope 
functions on the scale of the constituent bulk semiconduc- 
tor lattice constant. The zone-center SL wavefunctions 
and energies are found as described in Ref. la Sec. IIIB. 

Once the zone-center SL eigenstates and eigenener- 
gies have been determinted, the SL states for K 7^ 
can be obtained by application of a generalized 
SL K • p theory. A SL state |£,5,K) at finite K 
is expressed in terms of the zone-center SL states 
as (r|£,5,K) = exp(zK • r) ^^^^ C£5^s(K)(r|iV, 5, 0), 
where CcsNsCi^) are the expansion coefficients. By in- 
serting {r\C,S,K) into Eq. (O, multiphng the resulting 
equation on the left-hand side by exp(— K • r)(A^, S, 0|r) 
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and integrating over a SL unit cell, we have 
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where Pjvsw'S'W = {N, S,0\p\N' , S' ,0). 

Spin splitting in zincblende QW's can be described 
with an effective momentum-dependent internal mag- 
netic field, H(/:,K) = T-^/^fl{C,K), where T = 
(g/LtB?i~^)^, is the Bohr magneton, and g is the elec- 
tron g-factor. The components of the spin precession 
vector r2(£,K) are given by the following relationsiS.: 
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where the magnitude of the spin precession vector is given 
hyn = h-'\EcTiK)-EciiK)\. 

This spin splitting is due to the combined effects of 
spin-orbit interaction and spatial inversion asymmetry, 
where this spatial inversion asymmetry can arise from the 
bulk inversion asymmetry (BIA) of the constituent semi- 
conductors and the structural inversion asymmetry (SIA) 
of the QWs. The spin precession vector can be decom- 
posed into n{C,K) = 0(D)(£,K) + 0(i^)(£,K) accord- 
ing to its symmetry, where 11^^^(£,K) and fl^^'> {C,K.) 
are the momentum-dependent spin precession vectors 
due to bulk inversion asymmetry and structural inver- 
sion asymmetry, respectively. Although many analyses 
of these spin splitting fields focus on the K-linear terms, 
here we keep terms to all orders. 

In the linear approximations, both [df^\c, E)]"^ and 
[ft'"^\c, E)]'^ are assumed proportional to E for all en- 



ergies 
>(R) 



— . In addition, the nonlinear [rjg'^-' (C, E)] and 



[rig (£, i?)], are neglected in other analyses. Depend- 
ing on the size of the quantum well, these two terms can 
become dominant. The importance of these nonlinear 
effects was pointed out in symmetric quantum wells^^. 
Multiband calculations for asymmetric quantum wells 
have been done, but as a function of asymmetric dop- 
ing, not as a function of an applied electric fieldil^. 

Electron spin relaxation near room temperature is 
dominated by the D'yakonov-Perel' (DP) mechanism 
due to the effective internal magnetic field arises from 
spin-orbit interaction2iS*2£. In the motional narrowing 
regime the electronic spin system is subject to a time- 
dependent, randomly oriented effective internal magnetic 




" 0.0, 



50 100 150 200 

E (meV) 



50 100 150 200 

E (meV) 




- 1 




: ; 

L;i = 3j' 

: / 
: / 


/ n = 1 J 


: / 
- / 
: / 


/ L =150A: 
/ w : 




/ 1 lOOk^/cm: 



0.3. 



> 

CD 



0.2c>i 



0.1 



50 100 150 200 

E (meV) 



50 100 150 200 

E (meV) 



0.0 




200 400 600 

E (meV) 



200 400 600 

E (meV) 



FIG. 1: (color online). Energy dependence of electron 
spin precession vector for a 75-A and a 150-A GaAs/100- 

A Gao.6Alo.4As QW at 300K. [f^i;°'(l, E)]'^ as a function 
of E for a 75-A QW (a) and a 150-A QW (b) with 

T = kV/cm. [fl^^2il,E)f as a function of E for 
a 75-A QW (c) and 'a 150-A QW (d) with T = 100 
kV/cm. [d^;l{l,E)f as a function of E for a 75-A 
QW (e) and a 150-A QW (f) with T = 100 kV/cm. 



field H(£, E) which changes direction with an orbital 
scattering time t{C, E) that is much shorter than the 
precession time of cither the constant applied magnetic 
field Ho or the random field. The spin relaxation time 
and spin coherence time in zincblende QWs depend on 
the transverse H± {C, E) and longitudinal H\\ (£, E) com- 
ponents of the random field, according to 



dEV{£,E)Hl{£,E)T{C,E)h{C,E), {7a) 

r <J 



= Tf"E / dEV{C,E) 



-Hl{C,E)+Hl{C,E) 



T{C,E)h{C,E), 



(7b) 



where ^{C, E) is the density of states and h{C, E) is the 
electron distribution functioni^. 
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Now we present the results of our quantitative numer- 
ical calculations of the spin precession vectors, spin life- 
times, and spin coherence times of (001) GaAs/AlGaAs 
quantum wells. We find that the nonlinear effects are 
important and both [n'"r^\c, E)] and [r23^''(£, i;^)]are sig- 
nificant in these quantum wells. In the absence of an ex- 
ternal electric field, the effective internal magnetic field 
in a perfectly symmetric QW arises entirely from the BIA 
terms, and the SIA terms vanish identically. FiguresOJa) 
and n^b) show [rJx'^^(l, i?)]^ as a function of energy E 
calculated using Eqs. © for a 75- A and a 150- A QW, 
respectively, [fl'^i (l, E)]"^ is proportional to E and E^^, 
while [ri^^''(l, E)]"^ is proportional to E'^, consistent with 

the analytical resultsiiii^. [Q^n{l, E)]'^ decreases as 
the width of the QW increases; and it deviates from 
a linear dependence on E within 50 meV for QWs be- 
tween 50A and 150A indicating the importance of the 
nonlinear effects. 

When an external electric field is applied along the 
growth direction (^||z), the breaking of structural in- 
version symmetry results in modification of the effec- 
tive internal magnetic field through both the BIA and 
SIA terms [Figs. [Ti;c)-|Ti;f)]. [fl''^^ (1 , E)]^ increases as T 
increases, and it also increases with increasing [cf.. 
Figs. ^a)-^d)]. For the energy range of physical inter- 
ests, E)]"^ is approximately two order in magni- 
tude larger than [^'^^^{1, E)]"^; and E)]"^ deviates 
from a linear dependence on E for energies within 100 
meV [shown in Figs.^i^e) and^f)]. 

We find that the electric field dependence of r2^^]^(l, E) 
is linear, whereas that of fl''J^^{l,E) is nonlinear. The 
electric field dependence of fl^^-j E) is much weaker 

compared to that of Ql,^{l,E). To explore the electric 
field dependence of the effective internal magnetic field, 
ni^^(l,i?) and Q^n{i,E) as a function of T are calcu- 
lated for various values of E and the results for a 75-A 
and a 150-A QW are shown in Figs. |lla)-|2f)- For a 
constant value of E, the magnitude of fl''^^{l,E) and 

1 (1) £■) increases with increasing and the slope of 
the curves is positive for the BIA term, while it is neg- 
ative for the SIA term. The nonlinear effects are also 
manifested in the electric field dependence of E) 
[see Figs. |2i;a) and |a;b)]. 

Both Ti and T2 are more responsive to the electric 
field for wide QWs [i.e., [dTi^2/d!F) increases as in- 
creases], and this can be seen by examining the change 
in n^^n{l,E) with T [cf., Figs.[Ti;a)-|Ti;d)]. For example, 

the change in il^^j^'(l, E) from zero field to 100 kV/cm is 
much larger for the 150-A QW than for the 75-A QW. 
The calculations of electron spin lifetime and electron 
spin coherence time are performed using Eqs. iQ. To ex- 
amine the electric field dependence of spin lifetime and 
spin coherence time with applied magnetic field along the 
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FIG. 2: (color online). Electric field dependence of 
electron spin precession vector for a 75-A and a 150- 
A GaAs/lOO-A Gao.6Alo.4As QW at 300K. n''^}{l,E) 
as a function of T for a 75-A QW (a) and a 150-A QW 
(b). Vl''^1{1,E) as a function of T for a 75-A QW (c) 
and a 150-A QW (d). (e) as a function of T for a 
75-A QW. (f) as a function of T for a 75-A QW. 



growth direction, Ti and T2 as a function of T are calcu- 
lated for a 50-A, a 75-A, and a 150-A QWs at BOOK, and 
the results are shown in Figs.OIa) and|3Ib), respectively. 
For the DP mechanism, the spin lifetime and spin co- 
herence time are inversely proportional to the square of 
the effective internal magnetic field^"; and consequently 
both the spin lifetime and the spin coherence time de- 
crease with increasing electric field due to the increase 
of the total effective internal magnetic field with electric 
field. 

We find that the spin-lifetime enhancement originates 
from the destructive interference between the BIA and 
SIA effective internal magnetic fields due to the symme- 
try breaking of structural inversion symmetry in the pres- 
ence of an electric field. We also find that the calculated 
enhancement of the spin lifetime is much smaller than 
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FIG. 3: (color online). Electron spin lifetimes, co- 
herence times, and decoherence rates as a function 
of applied electric field for a 50-A, a 75-A, and 
a 150-A GaAs/100-AGao.6Alo.4As QWs at 300K with 
/i — 800 cm^/Vs assuming neutral impurity scatter- 
ing, (a) Ti as a function of with Ho||[001]. (b) 
T2 as a function of with Ho||[001]. (c) Ti as a 
function of J'^ with Ho||[110], Ho||[lTO], Ho||[100], and 
Ho 1 1 [010] for a 50-A QW. (d) Ti as a function of 
T with Ho 1 1 [110], Ho 1 1 [110], Ho 1 1(100], and Ho] [[010] 
for a 75-A QW. (e) Ti as a function of J- with 
Holl[110], Holl[lT0], Holl[lT0], and Holl[lT0] for a 150- 
A QW. (f) T^^ as a function of with Ho] [[001]. 



that of the previous calculations based on perturbative 
approachesi^, in which the zero-field spin lifetime is over- 
estimated and the nonlinear effects are not taken into ac- 



count. To investigate the electric- field-induced enhance- 
ment in the electron spin lifetime^^, Ti as a function of T 
for three QWs of different thickness is calculated for an 
applied magnetic field parallel to the in-plane direction 
and the results are shown in Figs.|3c)-|2Ic). As expected 
the calculated spin lifetime is identical for Ho|l[100] and 
Ho 11 [010] due to the crystal symmetry of (001) QWs, and 
Ti decreases with increasing J-. The electric field depen- 
dence of Ti for Ho 11 [110] is similar to that for Holl[100], 
and the spin lifetime is shorter for Holl[110] than for 
Ho 11 [100]. For Ho 11 [110], Ti increases as T increases, 
reaching a maximum value and then it decreases as J- 
further increases. The spin-lifetime enhancement factor 
[ri{T)/Ti{0)] decreases as increases, and the maxi- 
mum enhancement factor for a 50-A, a 75-A, and a 150-A 
QWs is approximately 3.8, 2.0, and 1.2, respectively. 

In the presence of an electric field, we find that spin 
decoherence depends on both the BIA and SIA effective 
internal magnetic fields, and the spin decoherence is dom- 
inated by the BIA terms at low field while it is dominated 
by the SIA terms at high field. To study the effects of 
the BIA and the SIA effective internal magnetic fields on 
spin decoherence (as opposed to spin lifetimes), the de- 
coherence rates (T^^) as a function of T for a 75-A QW 
due to BIA and SIA are calculated separately, and the 
results as well as the total decoherence rate are plotted 
in Fig.|3{f). It can be seen that the spin decoherence is 
dominated by the BIA effective internal magnetic field 
for an electric field as high as 100 kV/cm, and crossover 
occurs at approximately 150 kV/cm. The value of the 
crossover decreases as increases, and the crossover 
for a 50-A and a 150-A QWs occurs at approximately 
250 kV/cm and 120 kV/cm, respectively (not shown). 

In summary, we have studied quantitatively the elec- 
tric field dependence of electron spin coherence in (001) 
GaAs/AlGaAs QWs using a non-pertubative fourteen- 
band electronic structure theory. We find that both the 
BIA and SIA effective internal magnetic fields are electric 
field dependent. The electron spin lifetime and spin co- 
herence time are strongly influenced by an electric field. 
Even for moderate (< 100 kV/cm) electric fields, the spin 
decoherence is dominated by the BIA effective internal 
magnetic field. At optimal electric fields, the enhance- 
ment factor of spin lifetime for GaAs/AlGaAs QWs be- 
tween 50-A and 150-A is approximately 4, which is much 
smaller than that obtained from the previous calculations 
based on perturbative, linear approximations. 
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